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CALCULATIONS  OF  BRILLOUIN  BACKSCATTER 


IN  LASER-PRODUCED  PLASMAS 


I.  INTRODUCTION 

This  paper  studies  Briilouin  backscatter,  which  may  be  very  important  in  laser  fusion  schemes  in 
that  it  may  reduce  the  amount  of  laser  light  available  to  couple  to  the  pellet.  This  process  has  been  stu¬ 
died  experimentally  for  a  long  time.1'12  Experiments  in  which  the  plasma  has  a  short  scale  length  gen¬ 
erally  show  that  the  fractional  backscatter  saturates  as  the  irradiance  increases.  The  plasma  might  have 
a  short  scale  length  because  the  laser  pulse  is  very  short,  as  is  often  the  case  in  Nd  laser  produced  plas¬ 
mas.6  7  or  because  the  focal  spot  is  small  (say  less  than  20  wavelengths)  as  is  often  characteristic  of  C03 
laser  produced  plasmas.3  *  On  the  other  hand  long  scale  length  plasma  often  show  much  more  back¬ 
scatter  with  no  obvious  tendency  for  saturation  for  reflection  below  SO  or  60%.  One  way  to  produce  a 
long  scale  length  plasma  is  to  use  a  prepulse  to  produce  a  target  plasma  and  then  shoot  the  main  laser 
beam  into  it.  These  experiments  have  been  done  at  NRL8  9  recently  and  also  elsewhere,10  using  an  Nd 
laser  produced  plasma.  More  recently,  similar  results  were  found  with  a  C03  laser  produced  plasma." 
Long  scale  length  plasmas  can  also  be  produced  by  a  gas  jet.12  These  show  Briilouin  backscatter  which 
saturates  at  40  to  60%.  Often  in  laser  experiments,  the  fractional  backscatter  seems  to  be  almost 
uncorrelated  with  any  thing,6  7  that  is  graphs  of  fractional  backscatter  versus  irradiance  or  some  other 
parameter,  show  large  scatter.  One  possible  explanation13  is  that  very  small  prepulses,  which  might 
have  been  undetected  in  earlier  experiments  are  responsible  for  the  backscatter.  This  paper  is 
motivated  principally  by  experiments  in  which  the  backscatter  seems  to  depend  on  the  physical  parame¬ 
ters  in  a  fairly  smooth  and  predictable  way.  In  this  respect  the  NRL  prepulse  experiments  are  particu¬ 
larly  interesting  because  it  measures  the  density  profile  of  the  prepulse  plasma  up  to  a  density  of  about 
one-tenth  critical.*  9  Specifically,  it  finds  an  exponential  density  fall  off  with  a  scale  length  of  about  100 

Manuscript  submitted  March  '.  IdXI 

I 


MANHEIMER  AND  COLOMBANT 

fi.  The  density  profile  is  not  measured  above  0.1  but  the  scale  length  must  be  much  shorter 
there.13 

Our  method  of  calculating  the  backscatter  consists  of  numerically  solving  the  equations  for  the 
spatial  evolution  of  the  amplitude  of  the  laser  pump  wave,  reflected  wave  and  ion  acoustic  wave.  This 
approach  differs  from  other  similar  approaches14'16  in  that  we  do  not  use  a  plasma  slab  with  sharp  boun¬ 
daries.  Instead,  the  gradients  in  all  quantities  are  explicitly  accounted  for.  The  backscatter  ends  not  at 
any  artificial  boundary,  but  rather  when  the  plasma  density  gets  so  small  that  there  is  no  further 
amplification  of  the  reflected  wave.  Our  numerical  solutions  can  (and  often  do)  follow  the  interaction 
for  distances  of  hundreds  of  free  space  laser  wavelengths.  The  other  way  to  theoretically  study  Bril- 
louin  backscatter  is  particle  simulations.12  17  18  19  Recent  simulations18  have  in  fact  shown  some  aspects 
that  our  calculations  show.  The  problem  with  these  simulations  is  that  they  work  best  when  studying 
the  transient  response  to  the  backscatter  over  short  distances.  For  instance  a  simulation  of  1 04/  fl 
(about  10-  time  steps)  corresponds  to  only  5  pico-seconds  for  an  Nd  laser  produced  plasma.  Clearly  a 
simulation  with  length  and  time  scales  sufficiently  long  to  achieve  a  steady  state  for  a  plasma  with  real¬ 
istic  parameters  would  be  extremely  expensive.  Thus  the  type  of  calculations  we  do  here  appears  to  be 
a  necessary  link  between  studying  the  basic  physics  with  particle  simulations  and  interpreting  real  exper¬ 
iments,  especially  as  the  laser  pulse  length  and  plasma  size  increase. 

Section  II  sets  up  the  equations  for  mode  amplitude  and  presents  solutions  assuming  no  other 
nonlinear  effect  except  pump  depletion.  The  importance  of  velocity  gradient  and  electron  temperature 
in  determining  backscatter  is  demonstrated.  However  the  functional  dependence  of  R  (reflection)  on  / 
(irradiance)  is  very  sharp.  For  instance  a  factor  of  3  increase  in  /  can  increase  R  from  1%  to  50%. 

Section  III  examines  whether  nonlinear  effects  on  the  sound  wave  can  slow  down  the  dependence 
of  R  on  /.  There  are  two  nonlinear  effects  we  consider,  an  increase  in  the  sound  wave  damping,  and  a 
limit  on  sound  wave  amplitude.  As  in  Ref.  14,  we  find  that  an  increase  in  dissipation  has  some  effect 
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on  reflection  for  either  long  or  short  gradient  scale  length  plasmas.  However  this  effect  is  less  impor¬ 
tant  than  other  physical  effects  which  limit  reflection.  A  limitation  on  sound  wave  amplitude  however, 
does  have  some  effect  on  the  reflection  for  short  scale  length  plasmas,  but  does  not  for  long  scale 
length  plasmas.  Thus  a  nonlinear  limit  on  sound  wave  amplitude  is  thereby  identified  as  one  physical 
process  which  gives  rise  to  saturation  of  R  vs  /  curves  for  short  scale  length  experiments. 

Section  IV  presents  another,  much  more  economical  way  of  calculating  reflection  which  is  uni¬ 
formly  valid  if  the  sound  wave  is  strongly  damped,  and  gives  the  same  total  reflection  in  all  cases.  It 
can  be  incorporated  in  full  simulation  codes. 

Section  V  presents  calculations  for  a  wide  band  laser  pump  We  show  that,  as  discussed  by  oth¬ 
ers.  a  wideband  pump  can  significantly  reduce  reflection. 

Section  VI  presents  a  variety  of  calculations  of  Brillouin  backscatter.  Principally  it  again  addresses 
the  question  of  what  slows  down  the  increase  of  reflection  with  irradiance.  One  very  important  effect  is 
the  increase  of  electron  temperature  with  irradiance.  We  find  that  this  effect  alone  gives  rough  agree¬ 
ment  with  the  NRL  prepulse  backscatter  experiments.  A  combination  of  this  effect  and  a  nonlinear 
sound  wave  amplitude  limitation  appears  to  be  responsible  for  the  saturation  of  reflection  with  increas¬ 
ing  irradiance  for  short  pulse  experiments.  Also  Section  VI  presents  other  calculations  on  backscatter 
with  a  wide  band  pump,  with  subsonic  blowoff  velocity  and  with  spatially  oscillating  blowoff  velocity. 

Section  VII  suggests  several  experiments  which  can  test  the  concepts  developed  in  this  paper,  and 
finally  Section  VIII  presents  conclusions. 


II.  STIMULATED  BRILLOUIN  BACKSCATTER  IN  SUPERSONIC  FLOW 


The  equations  which  describe  stimulated  Brillouin  backscatter  in  a  spatially  varying  medium  are 


at 


-*■  V, 


y„  Ap  L’exp  i  f  <kr  -  kr 


kjiix 


(la) 
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Ar  -  y„  Ap  A'  exp  i  J*  (kp  -  kr  -  ks)  dx  (lb) 

jj  +  vp  +  vpj  Ap  -  ~yuAsAr  exp  /'  /  (kp-kr-k<)Jx  <lc) 

were  M2|  is  the  amplitude  of  the  action  density  of  the  mode  and  A  is  taken  to  have  the  phase  of  the 
electric  field  £.  That 

A  —  A  lx.  t )  exp  /  J*  kdx  -  (It  +  c.c. 

The  indices  p,  r,  and  s  denote  the  pump  (incident)  wave,  reflected  wave  and  sound  wave  at  frequencies 
ft,,,  H,  and  fts  and  the  k ’s  are  the  wave  numbers  of  the  associated  waves 


The  resonant  point  is  defined  as  that  point  in  the  plasma  k,  -  k%  +  kr.  For  Brillouin  backscatter 
ks  =  2 kp,  kr~-  kp.  If  the  flow  is  supersonic  toward  the  laser  cs  +  v  <  0  so  that  backscattered  light  is 
blue  shifted.  The  v ’s  are  the  linear  damping  rates  of  the  waves.  For  the  incident  and  reflected  light  we 
use  inverse  bremsstrahlung 

2 

^  np  Zln  A/3.5  x  105  T*i2 

while  for  the  sound  wave,  we  use  the  Landau  damping  rate 
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which  is  valid  for  ZT",,  »  7].  Notice  that  vs  is  a  rapidly  increasing  function  for  temp.,..uure.  The 
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Equations  (Ua,  b)  with  u  —  0  are  the  linear  equations  for  Brillouin  backscatter  discussed  by  Rosen- 
bluth.20  If  the  flow  is  supersonic  (that  is  v,  and  vs  have  the  same  sign)  these  equations  predict  spatial 
amplification  of  the  backscattered  wave  and  a  steady  state  treatment  is  quite  straightforward.  The  total 
amplification  4r  between  \  ”  — 00  and  \  =  +°°  is  given  by 

.•fr(°°) 


.4r(-°°) 


=  exp 


rryj  A r 
K  vr  vs 


(4) 


where  2 K  -  ~(k„  -  kr  -  ks).  If  the  flow  is  subsonic  (that  is  vr  and  v5  have  opposite  signs).  Eq.  (la. 
cbc  ' 


b)  predict  absolute  instability  in  uniform  bounded  media.  Rosenbluth  et  al.:i  have  shown  that  inhomo* 
geniety  changes  the  nature  of  the  instability  from  absolute  to  convective.  However  the  spatial  and  tem¬ 
poral  behavior  of  the  reflected  wave  is  quite  complex  it  is  not  apparent  from  Rosenbluth's  work  how 
one  could  do  a  steady  state  theory.  In  this  section  we  consider  only  supersonic  flow  and  deal  with  sub¬ 
sonic  flow  later.  Equation  (lc)  includes  the  effect  of  pump  depletion. 


The  basic  procedure  then  is  to  numerically  integrate  Eqs.  da.  b.  c)  from  somewhere  below  the 
critical  surface  toward  the  laser.  The  density  profile  is  as  measured  in  recent  experiments  on  Brillouin 
backscatter  with  structured  pulses. 
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with  /(  -  20  fx ,  !i  -  100  /a  and  .v„  the  location  where  =  1020  cm'3.  We  assume  here  a  Nd  laser  pro¬ 
duced  plasma.  These  results  can  all  be  scaled  to  other  laser  wavelengths  as  will  be  discussed  later.  The 
resonant  point  is  taken  at  -  0.12  n„  the  backscattered  irradiance  is  taken  as  10-8  times  the  incident 
irradiance.  Later  we  will  show  that  results  are  not  very  sensitive  to  these  choices.  Another  important 
parameter  to  choose  is  the  velocity  profile.  We  study  two  choices 

v(.v  )  =  -ac, 

(6) 

v(x)  =  -cs  1  +  — -  — 

n 

where  xtr  is  the  location  of  the  critical  density  and  a  is  a  parameter  greater  than  1.  The  former  case  has 
uniform  velocity,  the  latter  has  a  velocity  gradient  roughly  equivalent  to  that  of  an  isothermal  rarefrac- 
tion  wave.  For  all  cases.  ZTJ  T,  =»  10,  Z  =*  2.66  (CH?  target).  Plots  of  power  reflection  versus 
incident  irradiance  for  1  micron  laser  light  are  shown  in  Fig.  1  for  three  choices  of  electron  tempera¬ 
ture.  Tt,  -  1  keV,  3  keV  and  10  keV.  For  the  choice  of  parameters  (and  also  in  all  other  calculations 
we  do),  the  damping  of  the  incident  and  reflected  wave  is  negligible.  Also  shown  in  the  dotted  curve  in 
Fig.  1  is  the  functional  form  of  reflectance  versus  irradiance  as  predicted  by  Eq.  (7)  of  Ref.  14  for  the 
case  of  vs  (t',  in  the  notation  of  Ref.  14)  equal  to  a  constant.  It  is  clear  that  the  velocity  gradient, 
increases  by  about  an  order  of  magnitude  the  irradiance  required  for  a  given  reflection.  Notice  that  in 
all  cases,  reflection  increases  very  rapidly  with  irradiance  with  no  tendency  for  saturation  below  about 
70%  reflectance.  In  the  next  sections  we  discuss  physical  mechanisms  which  could  given  rise  to  satura¬ 
tion. 

III.  NONLINEAR  THEORIES 

There  has  recently  been  discussion  as  to  whether  nonlinear  effects  on  the  sound  wave  can 
significantly  limit  the  backscatter.  There  are  basically  two  effects  which  are  possible.  First  of  all.  there 
may  be  some  nonlinear  effect  which  limits  the  sound  wave  amplitude  4,;  and  secondly,  as  the  sound 
wave  grows  it  heats  the  ions  so  that  i/,  increases.  In  this  section  we  study  these  effects. 
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We  find  that  if  4S  is  limited,  the  backscatter  depends  crucially  on  what  its  phase  is  as  it  propagates 
away  from  the  interaction  region.  There  are  at  least  two  possibilities.  If  the  damping  and  coupling 
y„  are  weak,  it  oscillates  as  a  freely  propagating  sound  wave  so  that  -ls  ~  constant.  On  the  other  hand, 
if  the  damping  and  coupling  are  strong,  it  oscillates  as  a  driven  oscillator  As  —  exp/  J  (kp-  kr-  k$)  dx. 
For  the  former,  the  size  of  the  interaction  region  is  governed  by  the  phase  mismatch  scale  length  as 
discussed  by  Rosenbluth.:0  For  the  latter,  the  size  of  the  interaction  region  is  governed  by  the  dissipa¬ 
tion.22-3  The  amazing  thing  is  that  for  either  case  (in  the  linear  theory)  the  total  amplification  of  the 
backscattered  wave  is  as  given  in  Eq.  (4).  For  strong  damping,  there  is  a  weak  interaction  over  an 
extended  region;  for  weak  damping  there  s  a  very  localized  strong  interaction.  Both  give  the  same  total 
amplification. 
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Any  nonlinear  theory  which  relies  on  an  increase  in  vs  must  account  for  a  broadening  of  the 
resonant  region  also.  Of  course  the  problem  is  much  more  complicated  than  indicated  here  because  (a) 
as  increases,  the  resonant  region  ultimately  can  become  larger  than  the  plasma,  <  b>  pump  depletion 
occurs  and  (c)  variation  in  other  quantities  besides  wave  number  mismatch  must  also  be  accounted  for. 
Our  calculations  show  that  there  is  some  reduction  in  backscatter  resulting  from  an  increase  in  vs. 
However  this  effect  seems  less  than  the  effect  of  increasing  the  velocity  gradient  or  electron  tempera¬ 
ture. 


Next  we  consider  the  effect  of  a  limit  on  the  amplitude  of  the  second  wave.  Let  us  assume  that 
in  addition  to  an  increase  in  vs,  a  nonlinear  effect,  for  instance  ion  trapping,  actually  does  more,  and 
somehow  clamps  a  limit  on  the  sound  wave  amplitude.  Whether  this  is  an  important  effect  or  not  then 
depends  on  how  large  an  amplitude  the  sound  wave  would  achieve  if  there  were  no  limit.  The  shorter 
the  scale  length,  the  larger  the  sound  wave  amplitude  must  be  to  cause  a  given  reflection.  Thus  we 
expect  a  limit  on  sound  wave  amplitude  to  play  a  larger  role  for  short  scale  length  (i.e.,  short  pulse 
experiments)  than  for  long  scale  length  or  double  structure  (i.e..  long  or  double  structured  laser 
pulses). 

The  last  possibility  we  consider  is  that  the  sound  wave  amplitude  be  limited,  but  there  be  no 
increase  in  vs.  Thus  the  wave  propagates  like  a  freely  propagating  sound  wave.  Since  the  interaction 
now  weakens,  but  without  a  corresponding  increase  in  interaction  length,  the  backscatter  can  be 
severely  limited  in  this  case. 

We  examine  two  approaches  to  increases  vs.  The  first  is  just  a  straight  increase  of  this  quantity. 
The  solid  curves  show  in  Fig.  2  reflection  versus  irradiance  for  four  choices.  ZTJ  T,  *  10.  ZTJ  T.  =  6. 
v ,  ”  u> ,/ 2  and  v,  -  to,.  The  ratio  of  irradiance  needed  for  given  reflection  in  these  cases  are  roughly 
1 :2:3:4  so  that  ion  heating  can  give  rise  to  some  reduction  in  backscatter.  If  the  ion  heating  is  so 
intense  that  —  w,  (Z7j,  —  T,)  there  can  be  a  fairly  significant  reduction,  although  it  is  still  less  than 
either  the  effect  of  a  velocity  gradient  or  what  the  effect  would  be  for  a  homogeneous  slab  For 


t 
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instance  if  Ze  T/T,  =  10,  then  os/a>s  is  about  0.015,  In  a  homogeneous  slat,  an  increase  of  o,  from  this 
to  vjms  =  1  would  imply  that  for  the  same  leflectton,  the  irradiance  would  have  to  be  about  sixty  times 
greater.  However  for  the  inhomogeneous  plasma,  the  irradiance  is  only  about  4  times  greater,  so  the 
effect  of  an  increase  of  i>,  is  much  less  than  it  would  be  for  a  homogeneous  slab. 


Another  approach  to  the  problem  of  an  increase  in  is  to  postulate  some  nonlinear  dissipation  of 
the  sound  wave.  Thus  will  be  some  increasing  function  of  .-I,.  We  use  a  model  based  on  ion  trap¬ 
ping.  The  idea  here  is  that  as  the  wave  sound  grows  it  traps  ions  and  forms  a  non-ihermal  ion  tail  to 
the  distribution  function.  This  non-thermal  tail  Landau  damps  the  wave.  The  temperature  of  this  tail 
is  taken  to  be  the  ion  mass  times  the  square  of  the  wave  phase  velocity  or  roughly  ZT, ..  Once  the 
number  of  ions  in  the  tail  is  known,  the  Landau  damping  rate  can  be  calculated.  If  the  trapping  width 

is  Av  =  wnere  <b  is  the  amplitude  of  the  electrostatic  potential  fluctuation,  the  number  of 


ions  in  the  tail  is  J %  y;(v)dv  where  \n  is  the  wave  phase  velocity.  We  now  calculate  the  number  of 
ions  in  the  tail  of  a  Maxwellian  at  temperature  ZTC.  Doing  the  calculation  in  the  same  way  from  (v^  - 
Av)  to  oo,  the  number  of  ions  is 


n,  =*  n 


erf 


v,  +  Av 


( 2TJM ) 


1,2 


erf 


v„  -  Av 


(iTj.vn11 


erfc 


v,  ~  Av 


( 2ZTJM ) 


1/2 


(9) 


and  the  damping  rate  from  this  tail  would  be 


P.ST 


V_77_  m_  _ 

8  M  uH 


<D, 


1  +  z 


I.  2 


—  exp  —  1/2 


(10) 


0>f3(  1  +  bk;TJoi}Xt) 

with  oj,  defined  in  Eq.  (3c)  Thus,  as  Av  increases,  the  damping  rate  also  increases.  We  have  computed 
the  reflection  for  ZTJT,  —  10,  Tc  =  3  keV  using  Eq.  (10)  or  Eq.  (3c).  whichever  is  larger,  for  t\.  The 
results  are  shown  in  Fig.  (2).  the  dotted  curve  being  for  anomalous  v,  for  the  case  of  a  double  struc¬ 
tured  profile  with  a  velocity  gradient.  There  is  very  little  reduction  in  backscatter.  Although  the  back- 
scatter  is  not  significantly  reduced,  an  examination  of  the  detailed  structure  of  say  shows  that  the 
interaction  region  is  now  much  larger,  and  that  in  large  regions  of  the  i\,  A,  c,  <  1 
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Finally  let  us  note  that  this  change  in  the  character  of  the  interaction  as  i\  increase  has  been 
observed  also  in  particle  in  cell  simulations.18  When  an  ion  tail  is  generated  by  the  ion  acoustic  wave, 
the  strength  of  the  interaction  decreases  and  there  is  a  transition  from  what  Kindel  calls  coherent  to 
diffuse  scattering. 

The  work  of  Ref.  (14)  does  find  saturation  by  increasing  rs.  However  this  work  uses  a  slab  model 
with  sharp  boundaries.  The  width  of  the  slab  is  independent  of  »/s,  so  there  is  no  chance  for  the 
interaction  length  to  increase  as  v%  increases.  We  find,  on  the  contrary,  that  the  gradients  and  associ¬ 
ated  phase  mismatch,  play  a  very  important  role,  even  far  above  threshold. 

Let  us  now  consider  the  second  possibility,  namely  that  there  is  an  increase  in  dissipation  coupled 
with  a  maximum  value  of  —y  One  reason  for  such  a  limit  is  that  as  —y  tries  to  grow,  it  not  only 

Landau  damps  on  ions  in  the  tail,  but  it  must  also  drag  more  ions  into  the  tail.  This  could  limit  the 

amplitude  of  We  use  the  same  model  for  the  anomalous  dissipation,  but  now.  in  the  equations 

‘ e  J 

*  (?({) 

tor  and  .4,.  we  do  not  allow  -y  to  be  larger  than  0.05.  The  result  for  R  versus  /  is  shown  as  the 
dashed  curve  in  Fig.  2.  Note  that  the  reduction  in  reflection  is  very  slight.  The  reason  is  that  with  the 
strong  dissipation.  — -  dos  not  get  much  above  0.05  even  without  the  limit. 

*  C 

We  now  consider  the  third  possibility,  namely  that  as  the  ion  wave  traps  ion,  it  forms  a  coherent 
phase  space  vortex  structure  so  that  it  moves  as  a  freely  propagating  ion  wave  (perhaps  with  a  nonlinear 
frequency  shift2').  A  model  for  this  is  to  limit  the  ion  wave  to  some  maximum  amplitude  but  not  to 
increase  the  dissipation.  Thus  the  interaction  region  does  not  increase  and  if  it  is  small  to  begin  with, 
the  reflection  significantly  decreases.  In  Fig.  2.  the  dash-dot  curves  show  reflection  as  a  function  of 

^  (K 

irradiance  for  maximum  density  fluctuations  of  —y  =  0.05  and  -y  =0.1.  There  is  now  a  significant 
reduction  in  reflection  and  saturation  with  irradiance. 
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However  it  seems  unlikely  to  us  that  the  interaction  could  be  coherent  for  very  long.  The  ion 
wave  would  almost  surely  break,  and  the  trapped  ions  would  randomize,  leading  to  Landau  damping. 
This,  in  fact,  is  what  has  been  observed  in  recent  simulations.18  Also,  if  the  interaction  region  is  small, 
there  is  virtually  no  ray-retrace  in  the  backscattered  wave  since  it  requires  a  long  distance  to  set  up  the 
volume  hologram.26  27  Since  at  least  some  ray-retrace  is  a  feature  of  the  experiments,  a  freely  propagat¬ 
ing  ion  wave  with  coherently  trapped  ions  does  not  seem  to  be  a  viable  nonlinear  state. 


The  next  question  is  whether  by  shortening  the  plasma  length,  it  is  possible  to  achieve  saturation 
of  the  reflection.  We  have  examined  this  for  a  plasma  with  an  exponential  profile  with  a  density  gra¬ 
dient  scale  length  as  short  as  15  /a.  Using  the  dissipation  model  of  Eq.  (10).  we  find  only  very  slight 
change  in  the  shape  of  the  R  versus  /  curve  even  though  vj k%  c,  >  1  for  most  of  the  plasma.  Thus  a 

•  •  ^  (f) 

nonlinear  dissipation  alone  has  very  little  effect  on  the  backscatter.  However,  an  -jr  limit,  coupled 

with  an  anomalous  dissipation  does  have  a  significant  effect  on  the  backscatter  for  short  scale  length 
plasmas.  In  Fig.  3  we  show  R  versus  /  curves  for  a  30  ^  scale  length  plasma  with  a  velocity  profile 
given  by 


with 


•cs(.v  =  x„) 


1  + 


.V(T 


(11) 


/  -  30  n. 

The  solid  curve  is  the  result  for  no  anomalous  dissipation,  the  dashed  curve  is  for  anomalous  and 
the  dotted  curve  is  for  vf  given  by  Eq.  (10)  but  also  with  an  limit  of  0.05.  This  latter  curve  does 

‘  t' 

not  show  saturation,  but  does  show  a  much  slower  increase  of  R  with  /.  The  reason  is  that  in  the  short 

plasma  we  consider  here,  much  larger  values  of  —  than  in  the  longer  plasma  are  needed  to  produce 

*  «’ 

the  same  reflection.  In  Section  VI  we  will  see  that  if  the  variation  of  T,  with  /  is  accounted  for.  the 
increase  of  R  with  /  is  much  slower  still,  so  that  saturation  should  result. 
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IV.  CALCULATIONS  IN  SONIC  AND  SUBSONIC  FLOW 

Calculations  with  Eqs.  (1  a-c)  for  sonic  and  subsonic  flow  are  much  more  difficult  because  (a)  the 
equations  are  singular  at  the  sonic  point  and  (b),  for  purely  subsonic  flow  the  reflected  and  ion  acoustic 
waves  are  at  their  thermal  levels  at  different  points  in  the  plasma  so  one  cannot  easily  calculate  the 
reflection  wjih  a  single  integration  through  the  plasma.  However,  the  calculation  can  be  made  for  sonic 
and  subsonic  flow  if  one  uses  Eq.  (7)  for  -fs  and  then  derives  equations  like  Eq.  (8)  for  l.-Ll2  and  ! .-Ir | ’. 
The  resulting  equations  are  quite  straightforward  and  are 

V,  ±  M,|>  +  2U,P  %  +  W  A.!!  -  -  V,  UP  -  -  2v,  UP 


dx  '  cix 

-  2y,;  l.-Ll-V/vr 
( k „  -  kr  -  ks )2  +  (vjvs)2 


(12) 


Equation  (12)  are  quite  straightforward.  They  are  two  real  equations  and  can  be  integrated  from  the 
critical  surface  outward  to  calculate  the  reflection.  They  are  valid  everywhere  if  the  ion  wave  is  strongly 
damped.  If  the  ion  wave  is  weakly  damped,  they  give  the  same  total  amplification  as  Eqs.  ( 1  a-c)  for 
supersonic  flow. 


Notice  that  Eq.  (12)  depends  only  on  vs:  so  it  gives  the  same  amplification  for  positive  and  nega¬ 
tive  vs.  Thus  it  can  be  used  where  vr  and  vs  have  opposite  sign  (subsonic  flow)  and  it  gives  the  same 
amplification  as  found  in  Refs.  20  and  21.  Also  the  equation  is  no  longer  singular  at  v,  =»  0.  If  the 
velocity  is  subsonic  near  the  critical  surface  and  increases  uniformly  through  a  sonic  point,  so  that  it  is 
supersonic  at  lower  densities,  resonant  positions  with  >  0  (to  being  the  doppler  shifted  sound  wave 
frequency)  will  be  in  the  subsonic  plasma,  while  resonant  positions  with  fis  <  0  will  be  in  the  super¬ 
sonic  plasma.  If  co  -  0,  then  k ,  «  °°  at  the  sonic  point  so  that  a  sound  wave  in  the  supersonic  plasma 
cannot  propagate  into  the  subsonic  plasma  and  vice  versa.  Thus  for  a  plasma  with  a  sonic  transition, 
the  procedure  for  red  shifted  light  is  to  start  the  integration  of  Eq.  (12)  at  the  critical  surface  and  march 
outward  until  the  sonic  point  and  then  stop.  For  blue  shifted  light  the  procedure  is  to  start  at  the  sonic 
point  and  integrate  outward  into  the  supersonic  plasma. 
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We  have  compared  results  for  numerical  integration  of  Eq.  (12)  with  results  for  numerical 
integration  of  Eqs.  ( la-c)  for  both  strong  and  weak  damping.  In  all  cases,  graphs  of  reflection  versus 
irradiance  were  quite  close  to  each  other  and  had  basically  the  form  shown  in  Fig.  1.  However,  the 
detailed  spatial  structure  say  A ,  can  be  quite  different  depending  on  whether  the  damping  is  strong  or 
weak.  Finally  let  us  note  that  numerical  integration  of  Eq.  (12)  is  so  economical  that  it  can  easily  be 
incorporated  in  a  full  one  dimensional  fluid  simulation. 

V.  WIDE  BAND  INCIDENT  LASER  LIGHT 


The  equations  developed  in  the  previous  section  can  also  be  applied  to  steady  state,  nonlinear  cal¬ 
culations  if  the  laser  light  has  nonzero  frequency  bandwidth.  For  instance  say  there  is  a  single  reflected 
wave,  but  different  incident  waves  so  that  the  subscript  p  (on  say  ilp)  now  becomes  an  index  of  sum¬ 
mation.  Thus  the  equations  for  the  frequency  components  of  a  wide-band  incident  laser  source  are 


v,  ±  UP  +  2UP  ~  +  2vr  I.CP  -  -I  ';:2y\lAp  7: 

dx  dx  ~  (kp  -  kr  -  kt)2  +  iv J v%)2 


V,^I4P  +  2U,H^  +  2,,I.4.F 


2y,;U,N--l,;l-/v,: 
(kp~  kr~  k,)2+  (»>S/V,): 


(13) 


-•/t!  1  p\ 

The  reflected  wave  frequency  ilr  is  specified  and  there  are  now  different  pump  wave  frequencies  lir. 
The  sound  wave  frequency  for  each  pump  wave  is  given  by  Ur  -  O,.  Given  these  frequencies,  each 
wave  number  is  determined  by  Eqs.  <2a-c).  For  each  pump  wave,  there  will  be  a  resonant  location  in 
the  plasma  where  kp(.x )  -  krix)  +  k,ix). 


If  the  laser  irradiance  in  the  plasma  is  given,  then  X  l  -L-d})!2  is  determined.  However  if  there 

r 

is  a  broad  spread  in  frequencies  it  is  possible  that  only  a  very  small  portion  of  the  incident  spectrum 
resonant  with  any  reflected  wave  in  the  underdense  plasma.  Thus  the  backscatter  can  be  significantly 
reduced  with  a  broad  band  pump.  This  is  analogous  to  the  reduction  in  linear  growth  rate  of  the  Bril- 
louin  backscatter  instability  in  homogeneous  media,  by  a  broadband  pump.;s  As  a  rule  of  thumb,  one 
can  calculate  the  bandwidth  necessary  to  reduce  the  backscatter  in  an  experiment  as  follows.  Assume 
that  for  coherent  incident  light,  the  backscattered  light  has  bandwidth  5co.  If  this  8io  measures  the 
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spread  in  resonant  positions  for  the  scattering,  then  a  wide  band  pump  should  have  bandwidth  Au>  » 
6<u  to  significantly  reduce  the  backscatter. 

VI.  OTHER  CALCULATIONS  OF  BRILLOLIN  BACKSCATTER 

In  this  section  we  present  the  results  of  other  calculations  of  reflectivity  Since  most  graphs  of 
reflection  coefficient  versus  irradiance  have  the  same  form  as  those  shown  in  Figs.  1  and  2.  instead  of 
plotting  the  graphs  'hemselves.  we  simply  give  a  single  point  to  them 

The  graphs  themselves  can  be  easily  constructed  then.  For  most  of  our  calculations,  e  have 
used  Etts.  ( la-c).  However,  for  calculations  with  subsonic  flow  a  broadband  pump  and  flow  with  a  spa¬ 
tially  oscillating  velocity,  we  used  Eq.  (12).  Unless  otherwise  specified  ZT.  T  —  10.  the  density  profile 
is  as  given  in  Eq.  (4)  and  the  velocity  profile  is  either  as  given  in  Eqs.  (5a  or  b)  All  calculations 
assume  a  Nd  laser  produced  plasma.  Unless  otherwise  specified  —  0.12  ^  . 

A.  Dependence  on  Thermal  Source  and  Resonant  Position 

For  nearly  all  our  calculations  the  irradiance  of  the  reflected  wave  at  the  critical  density  is  smailer 
than  that  of  the  incident  wave  by  a  factor  of  10~8.  To  test  the  dependence  on  this  parameter,  we  have 
redone  the  calculation  taking  the  reflected  irradiance  smaller  by  10“*  for  the  case  where  the  velocity  is 
given  by  Eq.  (5).  We  find  that  the  reflection  coefficient  is  0.26%  at  an  irradiance  of  2.33  *  It)1'  W/cm: 
so  that  a  change  of  4  orders  of  magnitude  in  reflected  irradiance  changes  the  reflectivity  by  less  than  a 
factor  of  two  along  the  horizontal  axis. 

Next  we  have  examined  the  dependence  of  reflection  on  position  of  the  resonant  point.  We  have 
redone  the  calculations  with  resonant  positions  at  n  »  0.2  n 0.12  n  r.  0.08  and  with  velocity  gra¬ 
dient.  The  reflection  coefficients  are  respectively  10%  at  an  irradiance  of  1.35  x  10l?  W/cnr.  19%  at 
2.9  x  10I?  W/cnr  and  20%  at  3.5  x  101'  W/cm;.  In  addition,  we  have  calculated  the  x<  lection 
coefficient  where  all  four  reflected  waves  and  their  associated  sound  waves  are  simultaneously  in  the 
system  (but  where  the  initial  irradiance  of  each  reflected  wave  is  smaller  by  a  factor  of  4)  and  where 
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there  is  no  velocity  gradient  The  rellection  is  30%  at  an  irradiance  of  4  2  *  I014  W/cm:.  which  is 
nearly  the  same  as  for  the  single  wave  resonant  at  n.  -  0  12  //.,  and  for  no  velocity  gradient  Thus  the 
rellection  depends  very  weakly  on  both  the  resonant  position  and  also  on  whether  :here  is  a  single 
backscattered  wave  or  a  spectrum  of  waves  resonant  at  different  positions 

If  there  is  a  spectrum  of  backscattered  waves,  each  resonant  position  corresponds  to  a  different 
frequency  Thus  a  measurement  of  the  backscattered  spectrum,  especially  if  it  is  time  resolved,  can 
give  information  on  the  flow  velocities  in  the  blow-off  plasma 

B.  Reflection  Where  f  increases  with  Irradiance 

The  curves  in  Fig.  I  show  that  letlection  increases  very  rapidly  with  irradiance  ai  fixed  tempera¬ 
ture  This  rapid  increase  is  usually  not  observed  experimentally  One  explanation  is  that  the  tempera¬ 
ture  does  not  remain  fixed  as  the  irradiance  varies,  but  rather  increases  in  some  way  with  irradiance. 
We  consider  two  possible  scaling  laws.  T,  «  /'  \  and  T  «  /:  which  cover  those  generally  assumed 
Plots  of  reflection  versus  irradiance  for  the  two  scaling  laws  without  a  velocity  gradient  are  shown  in 
Fig  4  Figure  5  is  an  analogous  plot  with  a  velocity  gradient  In  the  former  case.  T  is  assumed  to  be 
1  keA  at  an  irradiance  of  8  <  10"  W/cnv  and  in  the  latter.  3  keV  at  2  *  10"  W'/cnv  If  T  <x  /;  , 
there  is  some  weakening  of  the  dependence  of  R  on  /  For  T  *  I1  ' ,  there  appears  to  be  a  saturation 
in  the  reflection  at  around  40  to  60  percent  at  least  over  the  range  in  irradiance  indicated  According  to 
linear  theory  ,  the  reflection  is  constant  if  //  T,  is  constant.  Thus  increasing  T,  can  greatly  reduce  the 
backscalter  at  given  irradiance 

Also  shown  on  Fig  5  are  measurements  of  reflection  versus  irradiance  taken  from  Fig  3  ol  Ref 
4  In  Ref  4.  the  density  profile  is  roughly  what  we  have  assumed.  Although  the  velocity  profile  is  not 
measured,  our  choice  of  an  isothermal  rarefaction  wave  is  certainlv  reasonable  We  recall  from  Section 
III  that  for  a  long  scale  length  experiment,  a  limit  on  sound  wave  amplitude  had  only  a  small  effect  on 
the  reflection  However,  the  agreement  between  experiment  and  theory  here  is  fairly  good,  so  that  the 


I  5 


M  ANHEIMER  AND  COLOMBANT 


increase  in  electron  temperature  with  irradiance  is  probably  responsible  for  slowing  down  the  increase 
of  R  with  /. 


For  short  pulse  experiments,  as  shown  in  Section  III.  a  limit  on  the  sound  wave  amplitude  does 
have  a  strong  effect  on  the  functional  dependence  of  reflection  on  irradiance.  It  seems  reasonable  then 
that  a  limit  on  sound  wave  amplitude,  coupled  with  an  increase  in  Te  with  /  can  give  rise  to  a  saturation 
in  reflection  versus  irradiance  in  short  pulse  experiments. 

C.  Calculations  with  Subsonic  and  Oscillating  Velocity 


Since  velocity  gradient  is  such  a  strongly  stabilizing  effect  on  stimulated  Brillouin  backscatter.  it 
seems  likely  that  one  could  reduce  the  backscatter  by  decreasing  the  velocity  gradient  scale  length.  One 
way  to  do  this  without  significantly  changing  the  flow-off  plasma  profile  is  to  superimpose  an  oscillating 
velocity  on  the  plasma.  This  could  be  achieved  by  oscillating  the  laser  intensity.  The  backscatter  can 
be  significantly  reduced  for  subsonic  flow  but  only  if  the  amplitude  and  wavelength  of  the  velocity 
modulation  are  carefully  chosen  It  the  wavelength  of  the  velocity  modulation  is  too  long,  there  is 
almost  no  effect  on  the  backscatter  Similarly,  if  the  wavelength  is  too  short  ihere  is  also  almost  no 
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effect;  presumably  because,  as  the  wave  amplifies,  it  averages  over  the  rapid  oscillations.  For  the  den¬ 
sity  profile  given  in  Eq.  (5).  we  find  the  optimum  wavelength  for  the  velocity  modulation  is  about  50  fi. 
Specifically,  choosing  the  velocity  profile 

v--0.3cs  1  +  —  +  v,  sin  —  x  (15) 

L  '  1 

with  -  107  we  find  a  reflection  of  10%  at  2.4  x  1015  W/cm:.  Thus  an  optimum  velocity  modulation 
can  increase  the  irradiance  for  given  reflection  by  nearly  an  order  of  magnitude.  However  the 
backscatter  for  the  oscillating  velocity  is  still  no  less  than  backscatter  for  the  supersonic  flow  profile  we 
have  examined. 

D.  Backscatter  for  Single  Structure  Profiles 


In  Fig.  <6)  are  shown  solid  and  dotted  curves.  They  show  the  irradiance  at  which  the  rellection  is 
50%  as  a  function  of  /  for  ZT{./'  T,  -  10  with  anJ  without  velocity  gradient.  VVe  have  also  plotted  as 
squares  the  irradiance  of  50%  reflection  for  /  -  30  jx  for  T,  -  6.6  keV  and  15  keV  a  50%  increase  in 
electron  temperature  has  a  fairly  large  effect  on  the  backscatter.  We  have  also  attempted  to  examine 
the  effect  of  an  oscillating  velocity.  However,  for  the  supersonic  flow  profiles  examined  here,  there  did 
not  appear  to  be  any  significant  reduction  in  backscatter 
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E.  A  Broadband  Pump 

Here  we  examine  Brillouin  baekscatter  for  broadband  pump.  The  density  and  velocity  gradients 
are  as  given  in  Eqs.  (4)  and  (5)  and  the  electron  temperature  is  T,  -  3  keV  and  ZTJT,  =■=  10.  The 
incident  spectrum  consists  of  5  or  9  waves  of  equal  amplitude  and  equal  frequency  spread.  The  center 
wave  always  has  a  resonant  position  at  n  -  0.12  ncr.  In  Fig.  7  are  shown  the  incident  total  irradiance 
for  10%  reflection  and  50%  reflection  as  a  function  of  bandwidth  All Here  AH,,  is  defined  as  the 
frequency  spread  between  the  center  frequency  and  maximum  (or  minimum!  frequency.  Note  that  as 
All,,  increases,  the  initial  tendency  is  for  the  baekscatter  to  slightly  increase.  The  reason,  apparently,  is 
that  as  the  resonant  positions  move  away  from  each  other  the  reflected  wave  is  resonant  with  the 
incident  spectrum  over  a  large  region,  but  not  as  strongly  with  any  wave,  in  such  a  way  that  the  net 
interaction  increases.  However  as  the  bandwidth  continues  to  increase,  the  reflection  begins  to  decrease 
quite  substantially.  The  reason  is  that  as  All  increases,  the  different  resonant  positions  move  out  of 
the  underdense  plasma  so  that  the  backscattered  wave  interacts  with  less  and  less  of  the  incident  spec¬ 
trum 

F.  Different  Laser  Wavelengths 

If  collisional  damping  of  the  incident  and  backscattered  wave  is  neglected,  all  equations  we  have 
used  are  invariant  to  the  scaling 

L  —aL 

v  —  V 

iu  — ■  wi  a  (18) 

n  —  nut * 

/  —  //a'. 

This  allows  one  to  scale  the  results  presented  here  to  different  laser  wavelengths.  However  when  scal¬ 
ing  to  shorter  wavelengths,  collisional  effects  become  more  and  more  important,  so  hat  the  reflection 
coefficients  should  be  progressively  smaller  than  predicted  by  the  scaling  law  in  Eq.  (18)  as  the 
wavelength  decreases. 
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VII.  SUGGESTIONS  FOR  FUTURE  EXPERIMENTS 

This  work  indicates  that  one  of  the  most  important  quantities  to  determine  in  the  blowoff  plasma 
is  the  velocity  and  velocity  gradient.  It  might  be  possible  to  measure  this  quantity  experimentally  by 
inserting  in  the  outer  edge  of  the  laser  target,  a  thin  layer  of  high  Z  material  which  would  line  radiate  at 
some  frequency.  This  transition  line  would  be  frequency  broadened  by  the  ion  temperature  and  fre¬ 
quency  shifted  by  the  plasma  flow.  The  shift  would  be  most  visible  if  looked  at  head-on.  but  the 
broadening  would  be  equally  visible  from  all  angles. 

When  the  laser  light  reached  this  material,  it  would  accelerate  outward.  The  idea  then  would  be 
to  get  a  time  resolved  measurement  of  the  Doppler  shift  of  a  transition  line  looking  head-on  from  the 
time  it  first  appears  until  it  disappears.  This  would  give  the  velocity  as  a  function  of  time,  and  integrat¬ 
ing,  would  give  the  velocity  as  a  function  of  space.  It  is  also  important  to  note  that  the  temporal 
motion  of  interferometry  fringes  does  not  give  information  about  the  flow  velocity,  but  rather  about  the 
motion  of  the  constant  density  surfaces.  The  two  are  essentially  unrelated  in  laser  produced  plasmas. 
Kruer14  and  others15  suggest  that  Brillouin  backscatter  gives  rise  to  strong  ion  heating  so  that  ultimately 
ZT„  ~  T,.  It  would  certainly  be  very  interesting  to  attempt  to  measure  the  ion  temperature  in  an 
underdense  plasma  in  which  Brillouin  backscatter  takes  place.  There  appear  to  be  two  possible 
approaches  to  this  problem,  both  of  which  could  indicate  ion  temperature  and  possibly  even  detect  the 
difference  between  heating  the  tail  and  body  of  the  ion  distribution  function.  The  first  is  a  variant  of 
the  method  we  have  just  discussed,  namely  seeding  an  outer  layer  of  target  with  a  high  Z  impurity,  only 
now  examining  the  Doppler  broadening  rather  than  Doppler  shift  of  the  associated  line  radiation. 

The  second  technique,  which  we  have  proposed  previously. ;s  is  to  seed  the  outer  part  of  the  target 
with  DT  or  D-He }  and  look  for  fusion  neutrons  or  protons  for  the  case  where  there  is  strong  Brillouin 
backscatter.  The  number  of  14  MeV  neutrons  or  protons  produced  depends  strongly  on  the  ion  tem¬ 
perature  and  Ref.  29  calculated  it  for  Maxwellian  ions.  It  should  also  be  possible  to  calculate  the  neu¬ 
trons  or  protons  produced  for  an  ion  distribution  with  a  non-thermal  tail.  Since  the  production  of 
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fusion  neutrons  is  most  likely  more  sensitive  to  non-thermal  tails  on  the  ion  distribution  function,  and 
the  Doppler  broadening  of  line  radiation  probably  most  sensitive  to  the  bulk  temperature,  a  comparison 
of  the  two  measurements  might  indicate  whether  the  ion  distribution  function  is  or  is  not  thermal. 

VIII.  CONCLUSIONS 

We  have  presented  here  a  large  number  of  calculations  of  Brillouin  backscatter  in  laser  produced 
plasmas.  We  find  that  at  given  irradiance,  the  backscatter  depends  quite  strongly  on  velocity,  velocity 
gradient  and  electron  temperature.  It  depends  more  weakly  on  ion  temperature  or  other  anomalous  dis¬ 
sipation.  For  a  short  scale  length  plasma,  the  backscatter  can  be  reduced  by  a  nonlinear  limit  on  the 
amplitude  of  the  sound  wave.  However  this  effect  becomes  less  and  less  important  as  the  scale  length 
increases.  The  saturation  of  backscatter  with  irradiance  characteristic  of  short  pulse  experiments  seems 

to  be  a  combination  of  saturation  by  -=r  limit  as  discussed  in  Section  III  and  increase  of  T,  wuh  irradi- 
ance  as  discussed  in  Section  VI.  The  lack  of  saturation  with  irradiance  in  double  pulse  experiments 

fc*  (i) 

results  from  the  fact  that  an  limit  no  longer  plays  an  important  role. 

'  t’ 

All  our  studies  were  done  for  Nd-laser  produced  plasma.  However  as  long  as  collisional  damping 
of  the  incident  and  reflected  light  is  not  important,  the  calculations  scale  easily  to  any  wavelength.  The 
two  principal  methods  available  to  reduce  the  backscatter  seem  to  be  in  control  of  the  blow-off  velocity 
(possibly  including  setting  up  a  spatially  oscillating  velocity),  and  increasing  the  bandwidth  of  the 
incident  light. 
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Table  1  —  The  Density  Gradient 
Scale  Length  at  the 
Assumed  Electron  Temperatures 
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R(%) 
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Kit  :  -  Double  siructured  density  protile,  solid  curves.  R  versus  /  Tor  £Tf  T  -  lit  and  n.  r,  w,  -  I  :  and  I.  doited  curve.  R 
versus  /  with  anomalous  dissipation:  dashed  curves  R  versus  /  for  anomalous  dissipation  and  ~  <  0  05:  dashed-dot  curves.  R 
versus  /  lor  —y  <  0.0?  or  0  !  but  with  no  dissipation 
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